Threefold Reduction in Thermal Conductivity of Vycor Glass Due to Adsorption of 

Liquid 4 He 
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We report thermal conductivity measurements of porous Vycor glass when it is empty and when 
the pores are filled with helium between 0.06 and 0.5 K. The filling of liquid 3 He and liquid 4 He inside 
the Vycor pores brings about respectively two and three fold reduction of the thermal conductivity 
as compared with empty Vycor. This dramatic reduction of thermal conductivity, not seen with 
solid 3 He and 4 He in the pores, is the consequence of hydrodynamic sound modes in liquid helium 
that greatly facilitate the quantum tunneling of the two-level systems (TLS) in Vycor and enhance 
the scattering of the thermal phonons in the silica network. 

PACS numbers: 67.25.de, 67.30.ef, 65.60.+a, 44.30.+V 
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The atoms or molecules in a disordered glassy mate- 
rial have multiple potential energy minima in their spa- 
tial configurations. At sufficiently low temperature, typ- 
ically below 1 K, the thermodynamic properties of the 
material are dominated by quantum tunneling between 
two adjacent accessible energy levels. The predictions of 
this two-level system (TLS) model QUI, 

including a spe- 
cific heat C scaling with temperature T and a thermal 
conductivity k scaling with T 2 , have been confirmed by 
experiments 

Zeller et al. found that the thermal conductivity of 
vitreous silica below 1 K is much lower than its crys- 
talline counterpart, quartz, and scales with T 2 instead 
of T 3 -Hi . Heat is conducted in solid by thermally acti- 
vated phonons and the thermal conductivity is given by 
k = (l/3)cphononul, where c p honon is the contribution to 
the specific heat per volume by the phonons, u the av- 
erage sound velocity and I the phonon mean free path. 
Normally c p honon~T 3 and u is weakly dependent on T. 
At low temperatures, phonons propagating in crystalline 
solid like quartz are primarily scattered at the defects and 
boundaries of the sample. Therefore I is temperature in- 
dependent and k^T 3 . For amorphous solid like vitreous 
silica, phonons are dominantly scattered by TLS that in- 
volve the coordinated motion and rotation of Si04 tetra- 
hedra groups. The physical dimension of each group, 
based on simulation and neutron scattering studies, is 
more than 1 nm TLS scatter phonons by ab- 

sorbing resonant phonons to tunnel to the excited state 
and re-emitting phonons to tunnel back to the ground 
state with a time constant r. The re-emitted or scat- 
tered phonons are generally propagating in a different 
direction from the incident phonons. TLS-phonon scat- 
tering results in an I that scales with T _1 |3|. As a result 
k^T 2 . If there are mechanisms to enhance the TLS tun- 
neling rate r _1 , more phonons will be scattered per unit 
time and k will be reduced further. We show in this pa- 
per that the infusion of liquid helium into Vycor pores 
dramatically reduces k. 

In addition to the TLS-phonon scattering, phonons are 



also scattered by the porous structure when propagat- 
ing in a porous glass like Vycor. Nitrogen adsorption 
isotherm measurements at T = 77 K show that Vycor 
has a porosity of 28% and a specific surface area of about 



100 m /g [Hi, Il2j. Transmission electron microscopy 
(TEM) revealed the pore space to be a network of multi- 
ply connected cylindrical channels [l2j ]. The diameter of 
the channel and the thickness of the silica "walls" , based 
on TEM and adsorption studies, are both found to be ^7 
nm [ll], Ell ■ The spatial correlation length of the silica 
structure (also the pore structure) in Vycor, £, is found to 
be distributed between 5 to 60 nm by small ang le X-ray 
and neutron scattering and TEM studies 12Hl4j . There- 
fore, phonons with wavelengths near or shorter than 60 
nm are expected to be strongly scattered by the Vycor 
porous structure because of the acoustic mismatch at 
the silica-pore interface. Phonons with wavelength much 
longer than 60 nm, on the other hand, will be insensi- 
tive to the porous structure and scattered solely by TLS. 
The spectrum of thermal phonons obeys the Planck's dis- 
tribution which shows a peak at Xdom = huj A.2ksT, the 
characteristic or the dominant wavelength [4j and the half 
maxima of the distribution at 0.5Acj om and SXdom- Here 
h is the Planck's constant, u the sound velocity and ks 
the Boltzmann constant. Since the temperature range of 
this experiment is between 0.06 and 0.5 K, the charac- 
teristic frequency of thermal phonons, fdom = 4.2fcsT'//i 
is between 5 and 40 GHz. Taking the transverse sound 
velocity in Vycor, uy veor to be 2200 m/s [lj|, Xdom is 
found to be 420, 170, 50 nm at T = 0.06,0.15,0.5 K, 
respectively. At low temperature when Xdom is much 
longer than £, phonons are mainly scattered by TLS. As 
a consequence, k should scale with T 2 . With increasing 
temperature, Xdom decreases towards £ and the scatter- 
ing effect by the porous structure becomes prominent. 
This will reduce n and result in sub-quadratic tempera- 
ture dependence. This behavior has been experimentally 
confirmed fTMl. 



In this Letter, we present the thermal conductivity of 
Vycor infused by solid helium, liquid helium and helium 



2 



heater 



Cu-Ni capillary 



diam. =3 mm 




T, 



copper base 
ixing chamber 



FIG. 1. Experimental setup of thermal conductivity measure- 
ments. 



films. Remarkably and counter-intuitively, a significant 
reduction in k, up to a factor of three, is observed when 
Vycor is infused with superfluid 4 He films, liquid 4 He and 
liquid 3 He. 

Fig. Q] shows our experimental setup. The Vycor rod 
is 3 mm in diameter and 22.8 mm long. It is secured me- 
chanically and thermally into a copper base with epoxy 
resin (Stycast 2850FT) and attached to the mixing cham- 
ber of a dilution refrigerator. Outside the copper base, 
the Vycor rod is sealed only by a thin layer of epoxy 
resin less than 0.07 mm thick painted on its outer sur- 
face. The thin epoxy layer, impregnating the pores on 
the surface, is able to hold pressure up to at least 85 
bar below 4 K but does not measurably contribute to 
the thermal conductance of the sample. Such an experi- 
mental configuration was recently used in studies of su- 
persolidity in solid helium 1|| 20j. Helium is introduced 
into Vycor through the copper base with a thin Cu-Ni 
capillary. A heater is attached at the top of the Vycor 
rod. A germanium thermometer is secured 12.7 mm away 
from the copper base reading T\ and another thermome- 
ter is attached directly onto the copper base reading T 2 . 
The copper base is at a uniform temperature because the 
thermal conductivity of copper is 10 times higher than 
Vycor 2l|. Thermal conductivity of the Vycor rod is 
given by K = LQ/AAT, where L is the distance between 



the germanium thermometer and the copper base, A is 
cross section area of Vycor, Q is the steady state heat 
power and AT = T\ — T%- Measurements are made by 
imposing a dc power between 0.1 and 1.6 nW to main- 
tain AT small but experimentally significant (typically 2 
mK). The measured k of empty Vycor from 0.06 to 0.5 
K is shown in both Fig. [2] and Fig. [3] As noted above, 
the deviation below the T 2 dependence for T > 0.15 K is 
clearly evident. Our results are consistent with previous 
experiments within a few percent 

After the measurements of empty Vycor, solid, liquid 
and adsorbed films of both helium isotopes were intro- 
duced into the Vycor pores. For ease of description, these 



samples will be identified as solid, liquid and film sam- 
ples. For clarity, we plot k of the solid and atomically 
thin film samples in Fig. [2] and k of liquid and super- 
fluid film samples in Fig. [3l When Vycor is infused with 
liquid or solid helium, the samples can be considered as 
composites consisting of the two intertwining (silica and 
helium) networks. However, the heat is still conducted 
primarily by the silica network. This is the case because 
the sound velocities (and hence phonon wavelengths) in 
liquid and solid helium are ~10 times smaller than that 
in silica. Therefore the phonon wavelengths in helium are 
always shorter than 60 nm for T > 0.06 K. As a result, 
phonon propagation along the helium network, if exists, 
is scattered by the porous structure much more strongly 
than that along the silica network. 

As shown in the inset of Fig. [3J K of the atomically 
thin 4 He and 3 He film samples are slightly (~3%) higher 
than that of the empty Vycor over the entire temper- 
ature range. The helium atoms in these thin films are 
immobile and tightly bound to the random silica surface 
with a typical adsorption potential higher than 25 K for 
the 4 He film (surface coverage — 13.6 /imol/m 2 ) and 
3 K for the 3 He film (n 3 = 22 ^mol/m 2 ) 22]. The 4 He 
film does not show superfluidity even at T — K. In our 
temperature range of interest, the adsorbed atoms are 
therefore inert and merely serve to thicken the silica net- 
work. These tightly bound atoms soften the acoustic mis- 
match at the silica-pore interface and lead to the slight 
enhancement of k as compared with the empty Vycor. 
A more noticeable enhancement of n is observed for the 
solid helium samples in the low temperature limit. How- 
ever, as T is increased above 0.2 and 0.36 K, k of the solid 
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FIG. 2. Thermal conductivity of empty Vycor and Vycor 
with solid 3 He and 4 He at 50 bar. k of solid 4 He sample at 
47 bar, not shown in the figure, was found to lie between the 
two data sets shown in the figure over the entire temperature 
range. Inset shows thermal conductivity of Vycor with inert 
helium films. 
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3 He and 4 He samples respectively crossover from being 
higher to being lower than the empty Vycor value. We 
attribute the reduction of k to the coupling between the 
phonons in solid helium and the TLS tunneling in silica. 
For empty Vycor, phonons absorbed by TLS in silica can 
only be re-emitted to the silica structure. But for Vy- 
cor filled with solid helium, TLS have a second channel 
to release energy by exciting phonons in solid helium via 
modulating the van der Waals interaction at the silica- 
helium interface [lH, 0] • The TLS- helium coupling fa- 
cilitates the TLS tunneling and enhances the tunneling 
rate r^ 1 . As a result TLS-phonon scattering in silica 
becomes stronger and n is reduced. However phonon ex- 
citation in solid helium is strongly suppressed in the low 
temperature limit when the phonon wavelength A# e ex- 
ceeds the pore dimension. This is why an enhancement 
of k is seen in the solid samples for T < 0.2 K due to 
the softening of acoustic mismatch across the silica-solid 
helium interface. With increasing temperature, Xue de- 
creases to become comparable and then shorter than the 
pore dimension thus opening the TLS-helium coupling 
channel and resulting in a reduction of k. The differ- 
ence in the crossover temperatures of the solid 3 He (0.2 
K) and solid 4 He (0.36 K) samples reflects the different 
transverse sound velocities Ut and hence Xue of the two 
solid helium samples. u t of 50 bar solid 3 He is 180 m/s. 
This translates to a A# e of ~10 nm to match the pore 
dimension of Vycor near 0.2 K. In comparison, u t of 50 
bar solid 4 He is 270 m/s which translates to a Xhb of ~10 
nm near 0.36 K. The modest difference in k (< 20%) be- 
tween the solid samples and the empty Vycor over the 
entire temperature range of the experiment confirms the 
reasoning stated above that heat conduction is primarily 
along the silica network and perturbed by the presence 
of helium. 

Although the density and (first) sound velocity of liq- 
uid helium are similar as those of solid helium, the mea- 
sured k of the liquid samples are dramatically different. 
Instead of a modest change, k of full pore liquid 4 He sam- 
ple is reduced threefold as compared with empty Vycor 
over the entire temperature range as shown in Fig. [3l A 
smaller but still sizable (a factor of 1.5) reduction of k 
is also seen in unsaturated superfluid He film samples 
(ri4 = 51.7, 65.3 ^imol/m 2 ). Interestingly, a two-fold 
reduction of k is found when the pores are filled with 
non-superfluid liquid 3 He. This suggests the dramatic 
reduction of k is a consequence of fluidity and not super- 
fluidity of liquid helium. 

The fluidity of the liquid helium enables hydrodynamic 
sound modes to be excited in the porous structure. Biot 
showed that such excitation requires the viscous pene- 
tration depth of the liquid S = \Jr)/irpf to be smaller 
than the pore radius D/2, so that a fraction of liquid can 
be decoupled and in relative motion with respect to the 
solid (in our case the silica) matrix [25|, [26| . Here r\ and p 
are viscosity and density of the liquid, / is the sound fre- 
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FIG. 3. Thermal conductivity of empty Vycor and Vycor with 
liquid 4 He, liquid 3 He and superfluid 4 He films with superfluid 
transition temperatures well above 1 K. 



quency. Similar to phonons in solid helium, the hydrody- 
namic sounds are also excited by TLS via the modulation 
of van der Waals interaction at the silica- helium interface. 
But in contrast to phonons in solid helium in the pores 
which must satisfy the boundary condition that their 
wavelengths should be comparable or shorter than the 
pore dimensions, there is no such wavelength (and hence 
frequency) restriction for the hydrodynamic sound modes 
because helium atoms are mobile. The sound modes con- 
form to and also extend to the entire multiply connected 
pore space. The frequency spectrum of this hydrody- 
namic sound is continuous and is given by the Planck's 
distribution of thermal phonons in silica with the addi- 
tional Biot condition: / > f c = 4r]/irpD 2 . Therefore the 
hydrodynamic sound is much more efficient in enhanc- 
ing t -1 of TLS in the silica over the entire temperature 
range as compared with phonons in solid helium. This 
is the key reason for the dramatic reduction of k in the 
liquid samples. The hydrodynamic sound is often named 
as slow wave or slow sound }27H3l| because its velocity is 
slower than the sound velocity in the solid matrix. John- 
son pointed out that the slow sound in superfluid 4 He 
confined in porous media is the well-known fourth sound 
[28| | . Fourth sound can be excited at any finite frequency 
within the Planck's distribution because the zero viscos- 
ity of superfluid gives f c = 0. In the case of unsaturated 
superfluid films where viscosity is also zero, third sound 
plays the role of slow sound in enhancing r _1 . The fourth 
and third sounds themselves do not contribute to thermal 
conductivity because superfluid carries no entropy. 

Slow sound can also be excited in liquid 3 He with a fi- 
nite viscosity as long as / > f c . The density of liquid 3 He 
at 3 bar (pressure of our sample) is 0.089 g/cm 3 . The vis- 
cosity ranges from 150 to 50 ^Poise between 0.15 K and 
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0.5 K 32j. As a result f c ranges from 4.4 to 1.5 GHz, 
substantially lower than the characteristic frequency of 
the slow sound /, equivalent to fdom, ranging from 12 
to 40 GHz in the same temperature range. Therefore, 
there is a significant fraction of liquid 3 He in the pores 
of Vycor that supports the slow sound to enhance t -1 . 
The slow sound in liquid 3 He is expected to vanish below 
0.07 K when the viscosity of liquid 3 He is large enough 
to completely lock the liquid to the silica matrix. Un- 
fortunately, the large heat capacity of 3 He and the low 
thermal conductivity of the Vycor prevent us from cool- 
ing the composite below 0.15 K. 

The TLS tunnelin g ra te in an amorphous system has 
been shown to be [3LI23I] 

M 2 E 3 , E 

■ COth -r^= (1) 



2TTh 4: pU 5 



2k B T 



He is larger than t s I_ by a factor of 4. As 
~ 1 is enhanced by a factor of 5. Since the 



Here E is TLS energy separation, p and u are the den- 
sity and sound velocity of the medium where TLS release 
energy, and M is the deformation potential that charac- 
terizes the coupling strength between TLS tunneling and 
the resultant strain in the medium. For empty Vycor, 
the deformation occurs at Si-0 bonds and Msi-o is on 
the order of 1 eV [|| . In this case the TLS tunneling rate 
is labeled as Tgl_ . With liquid helium in the pores, 
the TLS on the pore surface also couple with helium via 
modulating the silica-helium van der Waals interaction 
and Msio 2 -He is ^2 meV [23| . In this case the total 
TLS tunneling rate is given by r _1 = Tg^_ + Tg iQ _ He . 
Since the density and sound velocity of liquid helium are 
~10 times smaller than those of silica, Eq. ([I} shows 
that Tgl 02 _ 
a result, r 

physical dimension of the Si04 tetrahedra groups that 
make up the TLS is more than 1 nm and the thickness 
of silica "walls" is ^7 nm, at least 30% of the TLS re- 
side on the pore surface. The average r _1 of the entire 
sample is therefore enhanced at least by 2.2 times, close 
to the three-fold reduction in k found in full pore liquid 
4 He sample. This agreement supports the model we are 
proposing for the observed reduction of k. Supcrfluid 
films cause less reduction of n than full pore liquid 4 He 
because less helium are coupled to TLS. In the liquid 3 He 
sample, the slow sound propagates with finite damping 
because of the finite viscosity. The damping may dissi- 
pate energy to the silica network in the form of phonons. 
This may explain that liquid 3 He causes less reduction in 
k than liquid 4 He. 

Beamish et al. and Mulders et al. [IH also re- 

ported an enhancement of TLS tunneling rate in Vycor 
due to the adsorption of liquid and solid helium in ultra- 
sound experiments between 0.08 and 5 K. However, these 
experiments were not probing thermal phonons since the 
ultrasound frequencies correspond to characteristic tem- 
perature between 0.3 and 10 mK. Schubert et al. ob- 
served that the transmission of 25-GHz phonons from 



quartz to amorphous paraffin film is enhanced by the ad- 
sorption of liquid helium film on the paraffin at 1.8 K 
[33j. They attributed this to the enhancement of TLS 
tunneling in the amorphous paraffin film. 

In summary, we observe a dramatic reduction of ther- 
mal conductivity of Vycor when the pores are filled with 
liquid helium and superfiuid films. The reduction is a 
consequence of hydrodynamic slow sounds in liquid 4 He, 
superfiuid 4 He films and liquid 3 Hc that facilitate the 
TLS tunneling in silica and enhance the TLS-phonon 
scattering. 
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